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THE  CURRENT  CONVECTIVE  INSTABILITY 
AND  ITS  APPLICATION  TO 

DIFFUSE  AURORAL  SCINTILLATION  CAUSING  IRREGULARITIES 


INTRODUCTION 

Recently  data  from  the  DNA  Wideband  satellite  have  exhibited  high  latitude  (auroral  and 
subauroral)  scintillation  enhancements  which  appear  to  be  associated  with  north-south  gradients 
in  total  electron  content  [TEC;  see  Fremouw  et  ait  1977;  Rino  et  al .,  1978].  Rino  et  al. , [1978] 
have  shown  that  a regularly  occurring  scintillation  enhancement  can  be  identified  in  the  night- 
time auroral  zone  data,  in  the  region  of  diffuse  auroral  particle  precipitation  It  is  also  shown 
that  this  enhancement  is  due  to  sheet  like  F region  ionospheric  irregularities.  These  irregulari- 
ties occur  near  where  the  TEC  gradient  points  northward,  i.e.,  increasing  TEC  as  one  goes 
northward.  In  addition,  it  appears  that  a d.c.  electric  field  points  either  westward  or  northwest 
and  that  there  is  a very  shallow  plasma  density  gradient  in  altitude,  i.e.,  the  northward  TEC  gra- 
dient dominates  (C.  Rino,  private  communication,  1978).  The  diffuse  auroral  precipitation 
deposits  low  energy  particles  into  the  F region  without  any  E region  current  systems  being  set 
up. 

At  first  glance,  with  the  dominant  TEC  gradient  pointing  northward,  the  ambient  mag- 
netic field,  B0,  pointing  down  and  the  d.c  electric  field  horizontal,  this  would  appear  to  be  a 
prime  geometry  for  the  usual  F region  Ex  B gradient  drift  instability.  However,  with  the  d.c. 
electric  field  pointing  westward  or  northwest,  the  configuration  is  stable  (see  Fig.  1;  the  electric 
field  would  have  to  point  eastward  for  instability).  The  saving  feature,  however,  is  the  diffuse 
auroral  precipitation  (current)  which  acts  to  destabilize  the  above  geometrical  configuration  (see 
Fig.  1).  The  conditions  of  having  a current  along  B0,  a density  gradient,  Vn0,  perpendicular  to 
B0,  and  an  electric  field  perpendicular  to  B0  and  possibly  Vn0  can  result  in  instability  and  these 
types  of  instabilities  are  generically  called  current  convective  instabilities  [ Kadomtsev , 1965)  . 

This  type  of  instability  has  been  applied  to  the  positive  column  of  laboratory  gas 
discharges  [Lehnert,  1958;  Hoh  and  Lehnert,  1960;  Kadomtsev  and  Nedospasov,  1960;  Kadomtsev , 
1965],  but  to  our  knowledge  has  not  been  applied  to  the  conditions  in  the  F region  ionosphere. 
A simple  physical  picture  of  the  instability  is  as  follows  (see  Fig.  1).  First  let  us  discuss  the 
usual  F region  Ex  B gradient  drift  instability  picture.  In  the  upper  part  of  Figure  1 the  hor- 
izontal line  represents  an  unperturbed  contour  of  constant  electron  density,  the  background 
zero  order  electron  density  gradient  points  upward  (northward),  there  is  a horizontal  d.c.  elec- 
tric field  (eastward),  and  the  ambient  magnetic  field,  B0,  is  into  the  picture.  Now  let  the  den- 
sity be  perturbed  by  a small  amplitude  horizontal  sinusoidal  variation  (k^  parallel  to  E^.  The 
ions  Pedersen  drift  to  the  right  relative  to  electrons  whose  Pedersen  drift  is  for  all  purposes 
negligible  in  the  F region  (the  dashed  sinusoid  represents  the  electrons,  the  solid  the  ions). 
This  gives  rise  to  space  charges  (4-  and  — ) which  in  turn  cause  small  scale  electric  fields  E', 
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directed  alternately  to  the  left  and  right.  In  this  geometry  (eastward  E0),  the  corresponding 
E'x  B0  drifts  will  then  carry  the  enhanced  regions  downward  (southward)  and  the  depleted 
regions  upward  (northward)  so  that  they  both  will  appear  to  grow  in  amplitude  against  the  back- 
ground density  (convective  instability).  If  E0  were  to  point  westward,  as  in  the  diffuse  auroral 
case,  the  perturbations  would  disappear,  i.e.,  the  situation  would  become  stable  (also  if  Vn0 
were  reversed  in  Fig.  1 stability  would  result). 

Now  let  us  discuss  the  E x B geometry  with  a magnetic  field  aligned  current  as  illustrated 
in  the  lower  part  of  Figure  1.  There  is  an  ambient  electric  field,  E^,  perpendicular  to  B0  (point- 
ing down)  and  pointing  to  the  left  (westward),  the  density  gradient,  Vn0,  is  into  the  picture 
(northward),  the  current,  j_oH,  is  anti-parallel  to  B0,  the  horizontal  line  across  the  perturbation 
(where  k_  is  in  the  E^  — B0  plane  and  perpendicular  to  Vn0)  represents  the  unperturbed  con- 
stant electron  density  contour.  In  this  picture  without  j_o!,  the  system  would  be  stable  to  the 
Ex  B instability  (see  previous  paragraph).  The  projection  on  k of  the  ion  Pedersen  drift 
caused  by  E^  results  in  a drift  of  the  ions  along  k_  such  that  it  would  result  in  stability.  How- 
ever, the  direction  of  j_0n  implies  that  the  relative  drift  between  ions  and  electrons  (in  the  frame 
where  the  electrons  are  at  rest)  is  anti-parallel  to  B0.  This  motion  projected  on  k^  results  in 
establishing  space  charges  which  give  rise  to  small  scale  electric  fields  E"as  shown  in  Fig.  1 
(note  that  the  Pedersen  motion  due  to  E^  results  in  small  scale  electric  fields  opposed  to  E"). 
If  the  particle  motion  projected  on  k^  is  dominated  by  the  j_ol|  rather  than  E^  then  the  space 
charges  give  rise  to  total  small  scale  electric  fields,  E',  as  depicted  in  the  figure.  The  correspond- 
ing E'x  B0  drifts  will  then  carry  enhanced  regions  out  (southward)  of  the  figure  and  depleted 
regions  into  (northward)  the  figure  so  that  they  both  will  appear  to  grow  in  amplitude  against 
the  background  density.  From  this  picture  the  instability  criterion  is  |k  • Vd|  > |k_  ♦ (cEdi/B0) 
(i/j/flj) |,  where  J_ot|  — noeVd  and  vx  and  O,  are  the  ion-neutral  collision  frequency  and  ion 
gyrofrequency,  respectively.  This  will  be  mathematically  derived  in  the  next  section.  The 
lower  figure  can  be  made  similar  to  the  upper  figure  by  noting  that  j_on  = n0eVd  — o-,|E0n  and 
then  setting  ^0  Eq±  Hon.  We  note  that  in  the  above  picture  making  j_0u  in  the  direction  of  B0 
or  rotating _k  by  90°  about  E^  or  B0  results  in  stability.  We  also  note  that  the  above  picture  for 
instability  is  valid  even  for  E^  — 0.  If  E^  were  reversed  in  direction  the  parallel  current  would 
enhance  an  already  unstable  situation. 

This  instability  also  has  kM  « k±  (see  next  section)  and  so  the  irregularities  generated 
will  be  field  aligned.  This  current  convective  instability  can  directly  result  in  long  wavelength 
scintillation  causing  F region  ionospheric  irregularities  in  the  diffuse  auroral  region.  In  the  next 
section  we  present  the  theory  and  the  last  section  contains  the  summary. 

THEORY 

In  this  section  we  present  a theoretical  model  for  a long  wavelength  fluid  type  plasma  ins- 
tability which  may  account  for  the  scintillation  causing  diffuse  auroral  F region  ionospheric  irre- 
gularities. For  our  model  we  take  the  electron  density  gradient  to  be  pointing  northward  (y), 
the  ambient  electric  field  E0  is  in  the  westward  direction  (x)  and  the  magnetic  field  points 
downward  (z).  In  our  simple  model  we  have  equated  the  TEC  gradient  with  a gradient  in  den- 
sity. We  assume  that  the  horizontal  (northward)  electron  density  gradient  is  much  sharper  than 
the  altitude  density  gradient,  which  we  neglect.  Our  basic  equations  are 
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a“  + V • (n„V„)  - 0 

(1) 

c C c fi  i 

(2a) 

C C fie 

Ve  - -^-E  x z - f-  — 2.  E„  + Veoll 

B0  B0  vt 

(2b) 

V - j_  — 0,  j_  — InQqaVQ 

(3) 

where  the  subscript  a is  the  species  label  (e  is  electron,  i is  ion),  n is  density,  V is  velocity,  v is 
collision  frequency,  q is  charge,  E is  electric  field,  B0  is  the  ambient  magnetic  field,  z — B0/|B0|, 
Ha  - |qa|B0/mac,  j_  is  the  current  and  Va0  represents  the  diffuse  auroral  particle  precipitation 
velocity  along  B0  (M  and  ^ denote  parallel  and  perpendicular  to  B^  which  results  in  a zero  order 
current.  In  the  momentum  transfer  eqn.  (2),  we  have  neglected  inertial  and  temperature 
effects,  and  included  electron  collisions  to  first  order  for  completeness.  Equation  (2)  is  valid 
for  F region  ionosphere  altitudes  (vj£la  « 1).  In  (2a)  vx  is  taken  to  mean  ion-neutral  colli- 
sions, whereas  in  (2b)  vt  is  really  electron-ion  collisions.  We  have  neglected  the  electron 
Pedersen  drift  compared  with  the  ion  Pedersen  drift. 


We  assume  quasi-neutrality  so  that  ne  = n;  = n and  our  final  equations  are  the  electron 
continutiy  equation,  electron  and  ion  momentum  equations,  and  V • j_  — 0.  The  equations  are 
then  linearized  such  that  n -n0(y)  + n,E  - E0  x - V<£,  VQ  - Va0  + VQ  with  the  perturbed 
quantities  fi,  V ac  exp  i [kx  x *f  k|jZ  — cot],  where  to  = tor  + iy.  We  then  obtain 


^°ii  + ¥;iTE°*“tEoP 


Y eo  “ V eon  - -jj"  Eo? 


v s_  ik,^i  - -J-  -J-  ikx0x  + ikx0.p 

m^j  B0  1 1 j t>0 

V e — — - — ik,|0z  + ikx0y 


(4a) 

(4b) 

(5a) 

(5b) 


nle"e  B0 

Substituting  (4)  and  (5)  into  the  linearized  versions  of  (1)  and  (3)  results  in  the  following 
determinant  set  of  equations. 

i(w  - kl(Veol()  n -(■£-  ikx 


Bn 


8y  m.vc 


k,fno)0  - 0 


(6) 


t k*E°  i + Vdki 


fi  + nn 


mi*'i  meI/e 


+ -p  kx2  -£■ 

B0  fti 


0-0 


(7) 


where  Vd  =VioM  — VeoM  and  (6)  and  (7)  are  from  linearizing  the  electron  continuity  equation 
and  V • = 0 equation,  respectively. 
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From  (6)  and  (7)  one  obtains 


—l 


1 9^ 


^ By 


— E 

Bn  0 


n, 


+ vrf  — 


k l ne 

k*  "e 


C v ; 


B0Oi 


Eo  + — VH 


ili  + ik 

v,  Vr 


kii  v , 

kv2  n. 


where  w'  = w — k||VC0|,.  From  w = wr  + iy  we  then  obtain 


1 Bn0 


9y 


c c vi  , ,,  K|l 

b E«  7T7  + Va  k7 


£i  + ik 

Vi  V m 


n, 


(8) 


(9) 


which  for  the  lower  F region  » vJSlJ  becomes 


y 


J_  9rto 
n0  By 


C „ VjVc  K„ 

B0  0 nine  d kx 


^ c 

ae 


n,ne 


(10) 


We  see  that  in  (9)  or  (10)  y is  independent  of  |k|  and  only  depends  on  the  angle  that  k makes 
with  B c.  In  the  denominator  of  (9),  the  first  term  (in  brackets  ) multiplying  (kn/kx)T  comes 
from  the  parallel  motion  of  the  ions  and  electrons;  whereas,  the  remainder  of  the  denominator 
comes  from  the  ion  Pedersen  motion.  It  should  be  noted  that  the  instability  is  essentially 
unaffected  by  the  current  direction.  Thus,  downward  currents  work  iust  as  well  as  upward 
currents.  In  (9)  if  we  set  k|,  — 0 we  obtain  y = — (n  “^n^/Sy)  (cE^Bo)  which  is  the  usual 
result  for  the  Ex  B gradient  drift  instability.  For  our  geometry  this  shows  y is  negative  which 
implies  stability.  For  instability  in  (9)  y > 0 which  implies 


< 0 


(11) 


This  says  that  for  instability,  with  the  westward  E0,  we  must  have  Vd  kn/kx  < 0 and 
|Vd  kN|  > (k^Eo/Bo)  (i^/flj),  which  is  exactly  the  condition  set  forth  on  the  basis  of  our  physi- 
cal picture  presented  in  the  first  section. 


For  very  large  parallel  currents  (or  E0  — 0)  such  that  | Vdk M | » (v-JCl)  we 

have  from  (9) 


_L  9no  v kn 
n0  By  d kx 

^ W , 2 


n,  ne 

kii 

2 

V,  + ve 

k* 

(12) 


We  can  maximize  this  growth  with  respect  to  0 = kn/kx.  The  growth  rate  maximizes  for 
© = ± {y/niKn/z/,)  + ((l^e)]"1)1^  From  (9)  the  growth  rate  maximizes,  in  general,  for 
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0 = - (cE</B0Vd)  (vjil)  ± {(cEJB 0vd)2(^n,)2  + 

Wfls)  [(HA,)  + (n>e)]-r/2 


For  altitudes  — 350-400km  corresponding  to  the  observation  altitudes  [rtmo  et  al.,  1978]  vjily 
— 10  4,  vjElt  — 10'4  and  this  makes  kt/kx  = 0.7  x 10"4.  Consequently  for  maximum 
growth  the  growing  perturbations  will  be  nearly  perpendicular  to  the  magnetic  field.  From  (12) 
we  obtain  for  maximum  growth 


y max 


J_  !*°  V 

2n0  9y  d 


1 + 


-1/2 


(13) 


Equation  (13)  yields,  with  Vd 500  m/sec  (at  a density  n0  — 105  cm'3  this  yield  joM  — 8 ii 

amps/m2  which  is  in  line  with  experimental  measurements,  R.  Vondrak,  private  communica- 
tion, 1978)  L = noOnyOy)-1  — 50  km  (C.  Rino,  private  communication,  1979)  and  the  above 
collision  frequencies,  ymax  = 3.5  x 10"3sec"1.  For  a westward  electric  field  — 10  mV/m  (C. 
Rino,  private  communication,  1979)  the  growth  rate  is  diminished  somewhat.  For  this  case 
from  (9)  we  obtain  that  the  maximum  growth  rate  is  ymax  = 2.7  x 10-3  sec-1,  which  occurs  for 
k1(/kx  = 9.4  x 10"5.  Including  pressure  effects  in  the  problem  introduces  diffusive  damping  in 
(9).  A typical  cross-field  diffusion  coefficient,  is  — 0.2  m2/sec  and  a parallel  diffusion 
coefficient,  Dn,  is  — 108  m2/sec.  In  the  present  study,  these  effects  become  important  for  per- 
pendicular wavelengths,  Ai,  < 100  m and  parallel  wavelengths,  \M,  < 1000  km.  However,  typ- 
ical scintillation  causing  perpendicular  wavelengths  are  — 1 km  and  since  kn/kj^  — 10“4  we  are 
considering  highly  field  aligned  irregularities.  Larger  parallel  currents,  due  to  precipitation,  will 
of  course  produce  larger  growth  rates.  However,  too  large  a current,  Vd  > 1 km/sec,  would 
excite  the  collisional  electrostatic  ion  cyclotron  instability  [ Chaturvedi , 1976].  It  may  be  noted 
that  the  linear  theory  of  the  current  convective  instability  proposed  here  favors  a wavevector 
perpendicular  to  the  TEC  gradient  (as  well  as  1^). 


SUMMARY 


We  have  investigated  a simple  plasma  fluid  model  to  account  for  the  diffuse  auroral  scin- 
tillation causing  F region  ionospheric  irregularities  observed  by  the  DNA  Wideband  satellite 
[ Rino  et  al .,  1978],  By  taking  account  of  the  diffuse  auroral  particle  precipitation  (current)  the 
stable  Ex  B diffuse  auroral  geometry  (corresponding  to  the  observations)  becomes  destabilized 
by  this  parallel  current.  For  a westward  ambient  d,c.  electric  field,  E0,  and  a northward  dom- 
inant electron  density  gradient,  the  relative  drift  velocity  between  ions  and  electrons  parallel  to 
B0,  Vd,  must  satisfy  the  condition  — k * Vd  > k • (cE0/B0)  (^,/flj)  for  instability.  The  max- 
imum growth  rate  for  the  instability  is  y = n'^*  (9n0  / 9y)  Vd  [1  + (fle  i'i/fli*>c]“1/2/ 2.  The 
instability  is  mainly  field  aligned  (k,,  « kj).  The  instability  is  fluid-like  in  nature  and  so  can 
directly  account  for  the  long  wavelength  diffuse  auroral  scintillation  causing  F region  irregulari- 
ties. 
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Fig  I.  Simplified  physical  picture  of  E*B 
gradient  drift  instability  and  current  convective  instability. 
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BURLINGTON,  MA  01803 

OlCY  ATTN  RADIO  SCIENCES 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.  O.  BOX  983 
BERKELEY,  CA  94701 

OlCY  ATTN  J.  WORKMAN 


GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FDRGE  SPACE  CENTER 
GDDOARD  BLVD  KING  OF  PRUSSIA 
P.  O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OlCY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  D.  BDX  1122 
SYRACUSE,  NY  13201 

OlCY  ATTN  F.  RE  I BERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUOIES 
816  STATE  STREET  (P.D.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OlCY  ATTN  DAS I AC 
OlCY  ATTN  DON  CHANDLER 
OlCY  ATTN  TOM  BARRETT 
OlCY  ATTN  TIM  STEPHANS 
OlCY  ATTN  WARREN  S.  KNAPP 
OlCY  ATTN  WILLIAM  MCNAMARA 
OlCY  ATTN  B.  GAMBILL 
OlCY  ATTN  MACK  STANTON 
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general  electric  TECH  SERVICES  CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

01CY  ATTN  G.  MILLMfcN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  OIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  93111 

01CY  ATTN  JOHN  ISE  JR 
01CY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

Call  class  attn:  security  officer) 

01CY  ATTN  T.  N.  OAVIS  (UNCL  ONLY) 

01CY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

01CY  ATTN  TECHNICAL  LIBRARY 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  01 V 
77  A STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
URBANA,  1L  61803 

01CY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
URBANA,  IL  60680 

CALL  CORRES  ATTN  SECURITY  SUPERVISOR  FOR) 
01CY  ATTN  K.  YEH 


KAMAN  SCIENCES  CORP 
P.  0.  BOX  7463 
COLORADO  SPRINGS,  CO  80933 
01 CY  ATTN  T.  MEAGHER 

L INCAB  IT  CORP 
10453  ROSELLE 
SAN  01  EGO,  CA  92121 

01CY  ATTN  IRWIN  JACOBS 

LOWELL  RSCH  FOUNDATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL,  MA  01854 

01CY  ATTN  K.  BIBL 

M.I.T.  LINCOLN  LABORATORY 
P.  0.  BOX  73 
LEXINGTON,  MA  02173 

01 CY  ATTN  OAVIO  M.  TOWLE 
01CY  ATTN  P.  WALORON 
01CY  ATTN  L.  LOUGHLlN 
01CY  ATTN  0.  CLARK 

MARTIN  .MARIETTA  CORP 
ORLANOO  OIVISION 
P.  0.  BOX  5837 
ORLANDO,  FL  32805 

01CY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R . W.  HALPRIN 

oicy  attn  technical  library  services 


INSTITUTE  FOR  OEFEN5E  ANALYSES 
400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

OICY  ATTN  J.  M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLF HARO 
OICY  ATTN  JOEL  BENGSTON 

HSS,  INC. 

2 ALFREO  CIRCLE 
BEDFORO,  MA  01730 

OICY  ATTN  DONALO  HANSEN 

INTL  TEL  t TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  CAMINO  DEL  MAR 
OEL  MAR,  CA  92014 

OICY  ATTN  S.  R.  GOLDMAN 


MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


OICY 

ATTN 

P.  FISCHER 

OICY 

ATTN 

W.  F.  CREVIER 

OICY 

ATTN 

STEVEN  L.  GUTSCHE 

OICY 

ATTN 

0.  SAPPENF IELO 

OICY 

ATTN 

R.  BOGUSCH 

OICY 

ATTN 

R.  HENDRICK 

OICY 

ATTN 

RALPH  KILB 

OICY 

ATTN 

OA\E  SOWLE 

OICY 

ATTN 

F.  FAJEN 

OICY 

ATTN 

M.  SCHEIBE 

OICY 

ATTN 

CONRAD  L.  LONGMIRE 

OICY 

ATTN 

WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.  0.  BOX  208 
BEDFORD,  MA  01730 

OICY  ATTN  JOHN  MORGANS TERN 
OICY  ATTN  G.  HAROING 
OICY  ATTN  C.  E.  CALLAHAN 


/ 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JON  OASSOULAS 

LOCXHEEO  MISSILES  & SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  94088 

OICY  ATTN  OEPT  60-12 
OICY  ATTN  0.  R.  CHURCHILL 

LOOCHEEO  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  OEPT  52-10 
OICY  ATTN  RICHARO  G.  JOHNSON  OEPT  52-12 
OICY  ATTN  W.  L.  IMHOF  OEPT  52-12 


MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  OOLLY  MAO I SON  BLVD 
MCLEAN,  VA  22101 

OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELO  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  ATTN  E.  C.  FIELO  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 
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IONOSPHERIC  MOOELING  DISTRIBUTION  LIST 
UNCLASSIFIED  ONLY 


P 


PLEASE  DISTRIBUTE  ONE  COPY  TD  EACH  OF  THE  FOLLOWING  PEOPLE: 


ADVANCED  RESEARCH  PROoECTS  AGENCY  (ARPA) 
STRATEGIC  TECHNOLOGY  OFFICE 
ARLINGTON,  VIRGINIA 


HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MASS.  02138 


CAPT.  DONALD  M.  LEVINE 


DR.  M.  B.  MCELROY 
DR.  R.  LINDZEN 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

DR.  P.  MANGE 
DR.  R.  MEIER 

DR.  E.  SZUSZCZEWICZ  - CODE  7127 
DR.  TIMOTHY  COFFEY  - CODE  6700 
DR.  S.  DSSAKOW  - CODE  6780 
DR.  J.  GOOOMAN  - CODE  7560 


SCIENCE  APPLICATIONS,  INC. 
1250  PROSPECT  PLAZA 
LA  JOLLA,  CALIFORNIA  92037 


DR. 

D. 

A.  HAMLIN 

DR. 

L. 

LINSON 

DR. 

D. 

SACHS 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PENNSYLVANIA  16802 


DR. 

J. 

S.  NISBET 

DR. 

P. 

R.  ROHRS AUGH 

DR. 

D. 

E.  BARAN 

DR. 

L. 

A.  CARPENTER 

DR. 

M. 

LEE 

DR. 

R. 

DIVANY 

DR. 

P. 

BENNETT 

DR. 

E. 

KLEVANS 

UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES 

405  HILLGARD  AVENUE 

LDS  ANGELES,  CALIFORNIA  90024 

DR.  F.  V.  CDRONITI 
DR.  C.  KENNEL 


DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL  LABORATORY 
NOAA 

BOULDER,  COLORADO  80302 


UNI  VERS  I T>  OF  CALIFORNIA,  BERKELEY 
BERKELEY,  CALIFORNIA  94720 


DR.  A.  GLENN  JEAN 
DR.  G.  W.  ADAMS 
DR.  D.  N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 


DR.  M.  HUDSON 


UTAH  STATE  UNIVERSITY 
4TH  N.  AND  8TH  STREETS 
LOGAN,  UTAH  84322 


A.  F.  GEOPHYSICS  LABORATORY 
L.  G.  HANSOM  FIELD 
BEDFORD,  MASS.  01730 

DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.  M.  FORBES 
DR.  T.  J.  KENESHEA 
DR.  J.  AARONS 


OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 

DR.  H.  MULLANEY 


DR.  P.  M.  BANKS 
DR.  R.  HARRIS 
DR.  V.  PETERSON 
DR.  R.  MEGILL 
DR.  K.  BAKER 

CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK  14850 

DR.  W.  E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 
DR.  E.  DTT 

NASA 

GODOARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MARYLANO  20771 


COMMANDER  DR.  S.  CHANDRA 

NAVAL  ELECTRONICS  LABORATORY  CENTER  DR.  K.  MAEDO 

SAN  DIEGO,  CALIFORNIA  92152 


DR.  M.  BLEIWEISS 
DR.  1.  ROTHMJLLER 
DR.  V.  HILDEBRAND 
MR.  R.  ROSE 


PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NEW  JERSEY  08540 

DR.  F.  PERKINS 
DR.  E.  FRIEMAN 


U.  S.  ARMY  ABERDEEN  RESEARCH  AND  DEVELOPMENT  CENTER 


BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MARYLAND 

DR.  J.  HEIMERL 


INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY /NAVY  DRIVE 
ARLINGTON,  VIRGINIA  22202 

DR.  E.  BAUER 


COMMANDER 

NAVAL  AIR  SYSTEMS  COmAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 

DR.  T.  CZUBA 


UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA.  15213 

DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
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PHDTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  D.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.  d.  FREMOUW 

PHYSICAL  OYNAMICS  INC. 

P.  D.  BOX  1069 
BERKELEY,  CA  94701 

01CY  ATTN  A.  THOMPSON 

R & O ASSOCIATES 

P.  O.  BOX  9695 

MARINA  OEL  REY,  CA  90291 

0ICY  ATTN  FORREST  GILMORE 
0ICY  ATTN  BRYAN  GABBARD 
01CY  ATTN  WILLIAM  B.  WRIGHT  JR 
01CY  ATTN  ROBERT  F.  LELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  DRY 
OICY  ATTN  C.  MACOONALO 
OICY  ATTN  R.  TURCO 

RANG  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  EO  BEDftOZI AN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  ENO  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  O.  BOX  2351 
LA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  0.  SACHS 
OICY  ATTN  E.  A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUOBURY,  MA  01776 

OICY  ATTN  BARBARA  AOAMS 

SCIENCE  APPLICATIONS,  INC. 

HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

OICY  ATTN  DALE  H.  DIVIS 


SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 

OICY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  ORIVE 
PLEASANTON,  CA  94566 
OICY  ATTN  SZ 


SRI  INTERNATIONAL 
333  RAVENS WOOD  AVENUE 
MEM-0  PARK,  CA  94025 

OICY  ATTN  DONALD  NEILSON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  L.  L.  C06B 
OICY  ATTN  OAVIO  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  ATTN  RAY  L.  LEADA8RAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  0.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CDRP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  C SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

OICY  ATTN  R.  K.  PLEBUOH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  0.  OEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON,  MA  01803 

OICY  ATTN  CHARLES  HUMPHREY 
OICY  ATTN  J.  W.  CARPENTER 
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